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Figure 2. Examples are shown of Biomedical Image segmentation, morphing, and classification. 


GUI (see Figure 1) that allows users to 
do an end-to-end run of a patient’s 
record, including tooth segmentation of 
all teeth, grid morphing of each seg- 
mented tooth, and plaque classification 
of each tooth image. 

The automated and accurate process- 
ing of the captured images to segment 
each tooth [see Figure 2(a)] and then 
detect plaque on a tooth-by-tooth basis is 
a critical component of the PlaqTrak sys- 
tem to do clinical trials and analysis with 
minimal human intervention. These fea- 
tures offer distinct advantages over other 
competing systems that analyze groups of 
teeth or synthetic teeth. PlaqTrak divides 
each segmented tooth into eight regions 


using an advanced graphics morphing 
procedure [see results on a chipped 
tooth in Figure 2(b)], and a pattern 
recognition classifier is then used to lo- 
cate plaque [red regions in Figure 2(d)] 
and enamel regions. The morphing al- 
lows analysis within regions of teeth, 
thereby facilitating detailed statistical 
analysis such as the amount of plaque 
present on the biting surfaces on teeth. 

This software system is applicable to a 
host of biomedical applications, such as 
cell analysis and life detection, or robotic 
applications, such as product inspection 
or assembly of parts in space and industry. 

This work was done by Ashit Talukder, 
James Lambert, and Raymond Lam of Cal- 


tech for NASA’s Jet Propulsion Labora- 
tory. Further information is contained in a 
TSP (see page 1 ). 

In accordance with Public Law 96-51 7, the 
contractor has elected to retain title to this in- 
vention. Inquiries concerning rights for its com- 
mercial use should be addressed to: 

Intellectual Assets Office 

JPL 

Mail Stop 202-233 
4800 Oak Grove Drive 
Pasadena, CA 91109 
(818) 354-2240 
E-mail: ipgroup@jpl.nasa.gov 
Refer to NPO-3041 7, volume and num- 
ber of this NASA Tech Briefs issue, and 
the page number. 


@ Rapid Chemometric Filtering of Spectral Data 

Target species would be identified in real time. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


A method of rapid, programmable fil- 
tering of spectral transmittance, re- 
flectance, or fluorescence data to mea- 
sure the concentrations of chemical 
species has been proposed. By “pro- 
grammable” is meant that a variety of 
spectral analyses can readily be per- 
formed and modified in software, 
firmware, and/or electronic hardware, 
without need to change optical filters or 
other optical hardware of the associated 
spectrometers. The method is intended 
to enable real-time identification of sin- 
gle or multiple target chemical species 
in applications that involve high- 
throughput screening of multiple sam- 
ples. Examples of such applications in- 
clude (but are not limited to) 
combinatorial chemistry, flow cytometry, 
bead assays, testing drugs, remote sens- 
ing, and identification of targets. 

The basic concept of the proposed 
method is to perform real-time cross- 


correlations of a measured spectrum 
with one or more analytical function (s) 
of wavelength that could be, for exam- 
ple, the known spectra of target 
species. Assuming that measured spec- 
tral intensities are proportional to con- 
centrations of target species plus back- 
ground spectral intensities, then after 
subtraction of background levels, it 
should be possible to determine target- 
species concentrations from cross-cor- 
relation values. Of course, the problem 
of determining the concentrations is 
more complex when spectra of differ- 
ent species overlap, but the problem 
can be solved by use of multiple analyt- 
ical functions in combination with 
computational techniques that have 
been developed previously for analyses 
of this type. 

The method is applicable to the de- 
sign and operation of a spectrometer in 
which spectrally dispersed light is mea- 


sured by means of an active-pixel sensor 
(APS) array. The row or column dimen- 
sion of such an array is generally cho- 
sen to be aligned along the spectral-dis- 
persion dimension, so that each pixel 
intercepts light in a narrow spectral 
band centered on a wavelength that is a 
known function of the pixel position. 
The proposed method admits of two 
hardware implementations for comput- 
ing cross-correlations in real time. 

One hardware implementation would 
exploit programmable circuitry within 
each pixel of an APS array. The analog 
spectral-intensity reading of the pho- 
todetector in each pixel would be multi- 
plied by a gain proportional to value of 
the analytical function for the wave- 
length that corresponds to the pixel po- 
sition. As a result, the output from each 
pixel would be proportional to contri- 
bution of the pixel to the cross-correla- 
tion (plus background) . The outputs of 
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all the pixels along the spectral-disper- 
sion dimension would be summed to 
obtain the value of the cross-correlation 
(plus background) . 

Such on-chip cross-correlation could 
be performed rapidly because the ana- 
lytical function could be statically pro- 
grammed into the APS array and the 
multiplications could be done simulta- 
neously or nearly so. All of the addi- 
tions could be done simultaneously by 
means of a single binning instruction. 
The charge wells of all the pixels could 
be connected simultaneously, collect- 
ing all the charge outputs from multi- 
plication operations into one “super- 
pixel,” the single readout value of 
which would constitute the cross-corre- 
lation value for the given analytical 
function. For an instrument in which 
the APS rows were aligned along the 
spectral-dispersion dimension and in 
which the image of a spectrograph slit 
was aligned along the pixel columns 
and spanned multiple pixel rows, it 


would be possible to perform simulta- 
neous cross-correlations for multiple 
target species by applying, to each pixel 
row, the analytical function correspond- 
ing to one of the target species. A sepa- 
rate readout would be needed for each 
target species. 

In the other hardware implementa- 
tion, cross-correlations would be com- 
puted externally to the APS array. The 
multiplications and additions would be 
performed in pipeline fashion. If the 
APS-array outputs were analog, then 
programmable analog signals repre- 
senting the analytical functions would 
be synthesized in phase with the corre- 
sponding stream of analog APS-array 
outputs and the multiplications and ad- 
ditions would be performed by rela- 
tively inexpensive, commercially avail- 
able analog mixing and filtering 
circuits, respectively. If the APS-array 
outputs were digital, the cross-correla- 
tions could be computed by a digital 
signal processor. Ordinarily, the analog 


approach would be preferable because 
the analog operations can be per- 
formed much more rapidly than can 
the corresponding digital multiplica- 
tions and additions. 

This work was done by Gregory Bearman, 
Michael Pelletier, and Suresh Seshadri of 
Caltech for NASA’s Jet Propulsion Labo- 
ratory. Further information is contained in 
a TSP (see page 1 ). 

In accordance with Public Law 96-517, 
the contractor has elected to retain title to this 
invention. Inquiries concerning rights for its 
commercial use should be addressed to: 
Innovative Technology Assets Management 

JPL 

Mail Stop 202-233 
4800 Oak Grove Drive 
Pasadena, CA 91109-8099 
(818) 354-2240 
E-mail: iaoffice@jpl.nasa.gov 
Refer to NPO-30912, volume and num- 
ber of this NASA Tech Briefs issue, and 
the page number. 


Prioritizing Scientific Data for Transmission 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


A software system has been developed 
for prioritizing newly acquired geologi- 
cal data onboard a planetary rover. The 
system has been designed to enable effi- 
cient use of limited communication re- 
sources by transmitting the data likely to 
have the most scientific value. This soft- 
ware operates onboard a rover by ana- 
lyzing collected data, identifying poten- 
tial scientific targets, and then using 
that information to prioritize data for 
transmission to Earth. Currently, the sys- 
tem is focused on the analysis of ac- 
quired images, although the general 
techniques are applicable to a wide 
range of data modalities. Image prioriti- 
zation is performed using two main 


steps. In the first step, the software de- 
tects features of interest from each 
image. In its current application, the 
system is focused on visual properties of 
rocks. Thus, rocks are located in each 
image and rock properties, such as 
shape, texture, and albedo, are ex- 
tracted from the identified rocks. In the 
second step, the features extracted from 
a group of images are used to prioritize 
the images using three different meth- 
ods: (1) identification of key target sig- 
nature (finding specific rock features 
the scientist has identified as impor- 
tant), (2) novelty detection (finding 
rocks we haven’t seen before), and (3) 
representative rock sampling (finding 


the most average sample of each rock 
type). These methods use techniques 
such as K-means unsupervised cluster- 
ing and a discrimination-based kernel 
classifier to rank images based on their 
interest level. 

This program was written by Rebecca Cas- 
tano, Robert Anderson, TaraEstlin, Dennis De- 
Coste, Daniel Gaines, Dominic Mazzoni, Forest 
Fisher, and Michele Judd of Caltech for NASA’s 
Jet Propulsion Laboratory. Further informa- 
tion is contained in a TSP (seepage 1 ). 

This software is available for commercial li- 
censing. Please contact Don Hart of the Cali- 
fornia Institute of Technology at (818) 393- 
3425. Refer to NPO-40265. 


Determining Sizes of Particles in a Flow From DPIV Data 

The same equipment would be used to measure sizes as well as velocities. 

John H. Glenn Research Center, Cleveland, Ohio 


A proposed method of measuring the 
size of particles entrained in a flow of a 
liquid or gas would involve utilization of 
data from digital particle-image ve- 
locimetry (DPIV) of the flow. That is to 
say, with proper design and operation of 


a DPIV system, the DPIV data could be 
processed according to the proposed 
method to obtain particle sizes in addi- 
tion to particle velocities. As an addi- 
tional benefit, one could then compute 
the mass flux of the entrained particles 


from the particle sizes and velocities. 

As in DPIV as practiced heretofore, a 
pulsed laser beam would be formed into 
a thin sheet to illuminate a plane of in- 
terest in a flow field and the illuminated 
plane would be observed by means of a 
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